Van Hove Singularity and Superconductivity in Disordered 

Hubbard Model 



Grzegorz Litak 
Department of Mechanics, Technical University of Lublin 
Nadbystrzycka 36, PL-20-618 Lublin, Poland. 
(February 1, 2008) 

Abstract 

We apply the Coherent Potential Approximation (CPA) to a simple extended 
Hubbard model with a nearest and next nearest neighbour hopping for dis- 
ordered superconductors with s—, d- and p— wave pairing. We show how the 
Van Hove singularities in the electron density of states enhance the transition 
temperature T c for exotic superconductors in a clean and weakly disordered 
system. The Anderson theorem and pair-breaking effects in presence of Van 

Hove singularity caused by non-magnetic disorder are also discussed. 
Pacs. 74.62.Dh, 74.20-z 
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I. INTRODUCTION 



Magnetic and non-magnetic impurities in superconductors always attracted interest and 
their treatment played the essential role in theories of superconductivity. 



magnetic disorder on superconductors have been treated in many waysETcl. Their argu- 
ments, originally applied to classic BCS superconductors, were reexamined for novel, ex- 



Examining the influence of various kinds of magnetic and non-magnetic disorder caused 
by structural, substitutional, irradiation defects etc. on properties of conventional and un- 
conventional superconductors has shown that their responses to disorder are quite different. 
In contrast to conventional materials where only magnetic impurities affect the supercon- 
ducting properties, for unconventional superconductors the effect of both magnetic and 
non-magnetic disorder is usually strongH^S. 

Thus, it is not surprising, that the response to disorder become the fundamental criterion 
of unconventionality of the physical mechanism leading to superconductivity. Moreover, 
effects of disorder are of interest because the high temperature superconductors have to be 
doped (La 2 _a; Sr x Cu0 4 with Strontium x = 0.15, YBa 2 Cu30 7 __ a ; with Oxygen x = 0.1) to 
show the optimal critical temperature and doping is always accompanied by disorder. 

On the other hand in these quasi two dimensional layered materials the Fermi energy has 
been found in the vicinity of Van Hove singularities of the electron density of statesll~i3. This 
has lead to formulation of Van Hove scenario for high temperature superconductors which 
says that optimal critical temperature is reached when chemical potential passes through 
the Van Hove singularity in the density of state Doping with charge carriers does 

not only change the density in the system but also smear the density of states eliminating 
its singularities and, specially for anisotropic superconductors, introduce the electron pair- 
breaking phenomenon0fli@. 




non- 




s-wave, d-wave 
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Thus, one has to investigate the effect of rise and fall of the critical temperature T c near 
the Van Hove singularities very carefully, simultaneously, taking into account the effects of 
disorder. Clearly disorder reduces the critical temperature T c by elimination of singularities 
in the density of states and, at the same time, by braking pairs. The present paper is 
the extension of the previous ones0S@ where the electron hopping was introduced only 
between nearest neighbour lattice sites. Here following Refs. |68| , |69| we introduce the hopping 
to next nearest sites il (i.e. the ratio of hopping parameters t, t' for YBaCuO was suggested 
to be t'/t = 0.45) and investigate the effect of disorder on superconducting properties. 

In the normal states the nonzero hopping to the next neighbour sites t' introduces the 
distortion of the Fermi surface (Fig. la) which results in considerable changes in the electron 
density of states (Fig lb). Note, in case t' ^ the central Van Hove singularity produces 
the stronger enhancement of the density of states (N(E) for E « E vl ) than for t' = (Fig. 
lb). Moreover, the other Van Hove singularity associated with the bottom edge of the band 
E v2 creates the second peak in the density of states N(E). Thus, depending on the band 
filling, the additional electron hopping t' should have the effect on critical temperature T c , 
rising it to higher value. On the other hand, the distorted Fermi surface with the stronger 
dependence on k (Fig. la) seems to be less stable in the presence disorder as k is not a 
good quantum number in disordered system. 

The present paper is organized as follows. In Section II the discussion starts with the 
extended attractive Hubbard model defined on the square lattice which has the extended s- 
d— and p- wave order parameter solutions. Than we shortly overview and classify various 
types of Van Hove singularities in electron density of states of one band model and their 
influence on the superconducting critical temperature T c . In Section III we investigate dis- 
ordered superconducting systems and apply the Coherent Potential Approximation (CPA) 
for attractive Hubbard model. The Anderson theorem for non-magnetic impurity effect on 
isotropic s-wave solutions and the pair-breaking effect in case of anisotropic pairing is also 
discussed. Section IV is devoted to the examination of the disorder effect on the supercon- 
ducting critical temperature T c . Finally, Section V contains conclusions and remarks. 
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II. THE ROLE OF VAN HOVE SINGULARITIES IN THE CLEAN SYSTEM 



A. Bogoliubov de Gennes Equation 



We start with the single band Hubbard model with an attractive extended interaction 
which is described by the following Hamiltonian§!~@: 



H = E tiAaCja + \ E U i3 n i n i - 

ija ij i 



Si Hi. 



(1) 



In the above rii = n^+riii is the charge on site labeled i, fi is the chemical potential. Disorder 
is introduced into the problem by allowing the local site energy to vary randomly from 
site to site, c\ a and Cj CT are the Fermion creation and annihilation operators for an electron 
on site i with spin a, is the amplitude for hopping from site j to site i (with tu = 0) and 
finally is the attractive interaction (t/y < 0) which causes superconductivity and can be 
either local (i = j) or non-local (i ^ j). Starting from Eq. ([!]) we apply the Hartree-Fock- 
Gorkov§ approximation, which results in the Bogoliubov-de Gennes equation for a singlet 
(s- or d- wave) superconductor: 



E 



/ 
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£i — ^)^il — til — An 

-A* u (-Si + n)8u + tu 



^ (u \ 



E 
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(2) 



where and are electron and hole wave functions with up anti parallel spins j and j 
respectively The usual singlet one particle Green function, in the Nambu space G(i,j; iui n ), 
at the Matsubara frequency uj n = ^(2n + 1) satisfies: 



E 



lUJr, 



Si + n)8u + tu Aji 



^A*j {no n + Si - fj)Su - tu 



( 



Gn(l,j;zu n ) G 12 (l,j;iuj ri 



Sijl 



(3) 



/ 



G 21 (l,j;iw n ) G 22 (l,j;iuj v 

where the pairing potentials Ay will be taken to be non zero only when the sites % and 
j coincide (i = j), for on site interaction Ua, or are nearest neighbours, for off diagonal 
interaction Uij. On the other hand in case of triplet, p-wave, paring instead the Eq. ^| we 
have: 
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while the analogue of the equation of motion for the Green function in Eq. |3| should be 
rewritten as: 



E 



[{lU n -Ei + n)8n + tu] 1 



\ 



A+ 



[{iu n +Si- fi)8u - tu] 1 



G n (l,j;zuj n ) G 12 (l,j;iLu n ) 



Sijl 



(5) 



with the spin dependent 4x4 Green function. Each of its components G nm (l, j;iw n ) is 
defined as: 



G nm (l,j;zu u 



j;vjJ v ) 



n, m — 1, 2 . 



(6) 



The order parameter for triplet paring reads: 



A TT A Ti ^ 

A IT A II 



(r) 



We assume that the hopping integrals can take a nonzero values for the nearest and 
next nearest neighbours. For a clean system can be expressed in fc-space by Fourier 
transform: e k = J2 .t i:j e-' Ri i k as: 



= — 2t(cos fc x a + cos k y a) + 4t' cos k x a cos fc^a, 



(8) 



where t represents the nearest neighbour site amplitude of electron hopping, while t' corre- 
sponds to next nearest neighbour hopping, a denotes the lattice constant, ji is the chemical 
potential equal to Fermi energy in zero temperature (// = Ep for T = 0). We shall refer to 
the Greens function matrix as G(i,j;tu n ) which will be of 2x2 or 4x4 size for singlet or 
triplet solution. The above equations have to be completed by the self-consistency condition 
for pairing potential: 



= u^Yt^G^r^n), (9) 

P n 

for singlet pairing case and 

A f = UiAj2^ nVG i2(hr, ^n), a, a' =T, | (10) 

P n 

for triplet one, where rj is a positive infinitesimal, (3 = is the inverse of temperature T 
and Boltzman constant ks (in the units we use here fee = 1). To simplify matters we have 
assumed that the Hartree term Uy < nj- ff > can be absorbed into the hopping integral 
tij and dropped it from Eqs. and (UH). As usual Eqs. (§) and (|T0|) are to be solved 

together with the corresponding equations for the chemical potential fi that satisfy the 
following: 

n=lj2e^G u {i,i;iu n ) (11) 

P n 



for singles or 
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n = -Y i e«'nVGll(i,r,toj n ) (12) 



for triplets, where n is the number of electrons per unit cell. 

Here we do not wish to be very specific about the physical nature of the point defects 
represented by the site energies e*. We are rather going to provide a reliable analysis of the 
simplest possible non trivial model. Thus we take them to be independent random variables 
defined to have values \8 and — |d> with equal probability, |, on every site^B. As might 
be expected we shall be interested in the average of G(i,j;iu> n ) over the above ensemble. 
To calculate G(i,j; iuj n ) we shall make use of the Coherent Potential Approximation (CPA) 
which is the best method at hand for the mean field theory of disorder^!. 



B. Van Hove singularities in clean superconductors 

Let's assume that the sites form a square lattice. Then for clean system (£j = for all 
i), in the normal state, where Ay = 0, the spectrum is given by e& (Eq. |8|). It has a saddle 



point Van Hove singularity at E v = At', resulting in the logarithmic divergence of the density 
of states N(E) ~ — \n(E — E' v ) (Fig. lb) The density of states for a normal state N(E) is 
defined as: 

1 sr- ' ° 2 

where df is an element of Fermi surface (Fig. 2b) N(E) reaches its maximum if Fermi 
surface satisfy the relation E — e k = 0. Finally, density of states N(E) can be also expressed 
by the elliptic function of a first kind K(.E)S 



N ( E ) = lE^E-e k ) = f-J df\V k e k \-\ (13) 
iV u An z JE=e u 



ME) = — Y Im = 1 K 



16t 2 -(E- At') 



(14) 



iV^ E + zr]-e k 2 7rH^l + f^ \\ 16t 2 (l + f^) 

Van Hove singularities in density of states (Eq. 0,0) correspond to three characteristic 
flat regions of e k , where Vfcefc = 0. In Fig. 2a we have plotted band energy e k (Eq. ||). 
For cuprates the value of next nearest neighbour hopping term t' is usually considered as 
< t' < 0.5t. Here we have chosen t' = 0.45t (Eq. One can easily determine Van Hove 
singularities for saddle points: (|fc x |, |A;J) = (n/a,0), (0, ir/a) and these corresponding to 
band edges: the bottom (0, 0) as well as top one (it /a, it /a). Three isoenergetic contours '3', 
'2', T' have been marked in Fig. 2a for E/t = 2, 1.8, 1.6. They correspond to Fermi surfaces 
(Fig. 2b) for three values of band filling n = 0.338, 0.467, 0.584 respectively. Note that for 
n = 0.584, the Fermi surface has hole like characteristics while for n = 0.388, it corresponds 
to electron like system. For n = 0.467 the Fermi energy Ep = 1.8t passes through the Van 
Hove saddle point singularity. 

For the 'on-site' attraction (negative U) Ua = U the linearized gap equation (Eq. ||) at 
T c can be written as: 

i JL r dBi £ ^ -;>-*') tmh m\ = i r &e eei tanh m) , (15) 

W-oo iVV V 2 J 7T 7-oo 2E \ 2 J ' v ' 

where (3 C = l/(T c /cs), and T c is a critical temperature. 

If the interaction is off diagonal then the Fourier transform of Z7y leads to following 
expression: 



U(k — q) = —\U\ ( ^ k ^ q it ^ k ^Q _|_ 2 s i n k x a sin g^a + 2 sin /c y a sin g y a ) , ( I G ) 



where 



7 fc = 2(cos fc x a + cos k y a), 

r] k = 2(cos k x a — cos fc^a). (17) 

The pairing parameters for corresponding symmetry of solution: extended s- d- or p-wave 
one have the following form: 

Afc = Ao7 fc , 

Afc = Aor/ fc , (18) 
A£ = Aq sin(fc x a) + A v sm(k y a), 

and A^. is a Fourier transform of the matrix Ay (Eq. [7]). Despite of different types possible 
solutions described by Eqs. (BLBLO) the linearized gap equation for the critical temperature 
can be written in a compact form: 

"'"^ tanh (M) , (19) 



7T ./ x 2/-; \ 2 

depending on the solution symmetry a = s, d,p. 

The normal density of states N(E) = -k J2k ~ e fc) an d the corresponding projected 
densities (N S (E), N d (E) and N P (E)) used in Eq. (|T9|) can be expressed in terms of Green 
functions of the normal system as follows: 

N ( E ) = -^H -lmG u (k,E), 
N k n 



N k 4 71 

N d (E) = -±J2 ! jl^G 11 (k,E), (20) 
k 

N P (E) = -1^2(sinA; ;E a) 2 UmG xl {k,E). 
k 
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To calculate the above densities of states we have used the recursion method described 
in the Appendix A. The appropriate densities of states for the clean system N(E), N S (E) 
Nd(E) are plotted in Fig. 3a. Figure 3b shows the p-wave projected density N P (E) in 
comparison to N(E). 

The full line, in Fig. 3a, corresponds to the electron density of states N(E) with the 
spectral dependence of efc (Eq. 4) and t' = 0.45. The position of the central Van Hove 
singularity E v = — 1.8t (Fig. 2) corresponds to band filling n = 0.467. Apart from that 
singularity one can see another sharp peak in the lower edge of band and much smoother 
one in the upper edge. Both of them are singularities of band edges: bottom and top 
respectively. For the 'on site interaction' Ua the shape of N(E) in Fig. 3a and the gap 
equation (Eq. [L5D indicates that the critical temperature T c should be enhanced effectively 
around first two singularities for relatively low electron densities n < 1. As a matter of fact 
such situation can be seen in calculations of T c (Fig. 4, clean system for curves denoted 
by T'). The critical temperature T c obtained from the Eqs. ( p~9| ; [20| ) for various symmetries 



of the order parameter are depicted in Fig. 4a-c by lines denoted by '1' as a function of 
band filling n. Figure 4a corresponds to extended s— , d- wave cases while Fig. 4b shows T c 
to the p-wave solution. In this case there exists a maximum in projected density of states 
N P (E) close to Van Hove singularities in N(E) but it is not so sharp as in N(E) or N S (E) 
and Nd(E) because of the additional smearing term 2(sin k x a) 2 in the formula for N p (E) 
(Eq. ). Nevertheless, this secondary peak also produces the enhancement of critical 
temperature T c (Fig. 4b). For comparison, in Fig. 4c we present T c for an isotropic 'on-site' 
s- wave solution. Clearly, these calculations support the Van Hove singularity scenario. 
Note that for relatively low temperature Van Hove singularity is passing at n = 0.467 (Figs. 
2, 3a). The lack of the particle- hole symmetry results in the shift of maximum value of T c 
to higher value of electron densities so in the optimal doping the superconductor is hole 
superconductor. 

Moreover for n — > 2 we observe the degradation of critical temperature T c . In in this 
limit of slowly changing density of sates we can apply the result for the constant density of 
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states: 



N(E) 



1 



Q(E-D/2) Q(-E + D/2) 



(21) 



D 



where D denotes the bandwidth and 0(E) is the Heaviside step function. 

For small 'on site' attractive interaction U, the critical temperature T* is given by the 
analytic formula@0 



In our case, due to the asymmetry of the electron density of states N(Ep) << 1/D for 
n — > 2 the degradation of T c is faster T c « T*. 

The Van Hove singularities in the spectrum e k show up also in the projected densities of 
states N S (E), Nd(E) (Fig. 3a). Due to the factors 7 fc and rj k (Eqs. |l7l , [2"0| ) the dependences 
T c (n) are severely modified (Fig. 4a, curve '1'). Moreover, the maximum value of N p (E) 
(Fig. 3b) is also close to the Van Hove singularity and it results in the optimum of T c (Fig. 
4b, curve '1'). 

It is worthwhile to notice here that singularities in the band edges are important for an 
extended s-wave case while the saddle point one, located in the middle of the band, for 
d-wave pairing, similarly to an isotropic s-wave case. 

The positions of the Van Hove singularities result in the strong band filling n dependence 
of T c (Fig. 2b). Again the pairing dominates in selected regions of n. In case of d-wave it 
is the middle region of n ~ 0.5 while for extended s-wave pairing it is rather low electron 
densities region (n — > 0). The high electron density ( n — > 2) is also possible however the T c 
is much smaller. Same relation, as in the 'on-site' pairing case, governs the basic behaviour 
(Eq. |2^). However N(E) = 1/D should be substituted by the corresponding projected 
density of states N a (Ep). Thus, away the Van Hove singularity, (n — > 2): 



In all cases the Van Hove singularities play the mayor role and could be identified as the 
source of a critical temperature T c raise. Its dependence on doping n should be described 




(22) 




(23) 
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rather by strongly changing function contrasting with the case of a constant density of states 



However it is also worthwhile to note that depending on the paring symmetry different 
Van Hove singularities can matter and because of non symmetric character of densities of 
states this effect is approximate. Due to the this asymmetry we observe the additional 
shift of the optimal doping toward the center of a band. Interestingly, different Van Hove 
singularities and corresponding shifts depend on solution type and the interaction range. 
Eventually, an extend ed s-wave solution appears to be an electron superconductor while 
d-wave can be identified as hole one. 

The formulae for a critical temperature T c (Eqs. [19]) are based on the integral over the 
appropriate DOS which posses Van Hove singularities. The effect of shift is stronger for 
smaller T c where hyperbolic tangent is smearing the function under integrals. It is also 
worth to note that Van Hove scenario is working better for superconductors with relatively 
small transition temperature T c (which corresponds to a small interaction parameter U). 
This can easily be seen from the following function: 



present in the gap equations (|5|,[P2|). It can be interpreted as leading to a natural cut-off Eq 
around the chemical potential /i. Note that if the temperature T c is small then the function 
F(T C , E — fi) is non zero in the narrow range of energies around fi only. In fact in the limit 
T c — > it tends to the Dirac delta function (F(T C ,E — /i) — > 5(E — fi)) and the cut-off 
is limited to the neighbourhood of the E = fi point. Note, that for finite T c , Eq « 2T c . 
Consequently, for the logarithmic Van Hove singularity in the density of states near fi has 
the form: 



(Eq. where T c is simply Jn{2 — n). 





(24) 



N(E) « -N In 



E-fi 



(25) 



D 



and we get the Labbe-Bok formula for TcHHil: 




(26) 
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III. CPA FOR DISORDERED HUBBARD MODEL 



As mentioned above, the technical question we shall answer in this paper is what happens 
to above behaviour when the site energies £j are not the same on all sites but are randomly 
distributed. For example in a binary alloy, A c Bi^ c , we have random distribution of site en- 
ergies: £j = (ea, sb) depending on occupation at the site i by atom A or B with probabilities 
Pa,(b) — c and (1 — c) respectively. Such problem have been dealt with within the Coherent 
Potential Approximation (CPA) on a number of occasions in the pasl&^'El^iiS. Here we 
shall follow the usual arguments generalized as appropriate. In short, we shall take the CPA 
to mean that the coherent potential = in a site approximation, is defined 

by the zero value of an averaged t-matrix T(i, i; E). Namely 

< T a (i, i; E) > = £ P a T a (t, %■ E) (27) 

i 

[V a - £ CT (£))(1 - [V a - X(E)]G(i, i; E))- 1 ) = 0, 



where a = A, B specifies the occupation of the site i and hence the disordered potential V a . 

In case of 'on site attraction' (Eq. |l]with lly = Uudij). In the limit of small fluctuations 
in the paring potential = Aj), a constant averaged value 

Aj -> A for all i (28) 

can be applied. 

Using now CPA equations equations (Eqs. ^7| and Appendix B) it can be readily shown 
that due to the disorder A and E in the clean limit is renormalized to A and E in the same 
way 

~ x 2E - Tr£(£0) / s 

A{E) = A IE (29) 

- _ 2E - TrS(E)) 

2E 

On account of the general symmetry between the the averaged Green function G(i,i;E) 
elements and that of the coherent potential "S(E), for complex energies, it follows that 
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G u (i,j;E) = G 22 (i,j;-E*), 



(30) 



£ n (£) = £ 22 (-£*). (31) 



Note that one can write Eq. ( ^5| ) in terms of Matsubara frequencies (i? = zu; r , 



A(a) n ) uj n 



A w 



(32) 



n 



where 



iuj n = iuj n - £n(«a> n ) (33) 

Eventually Eq. ( pi]) leads to the linearized gap equation (Appendix B): 

/•oo tanh(^) , s 

l = Uj_JEN{E) (34) 

This equation relates directly to the similar one of the clean case (Eq. [15|) with the difference 
due to the substitution of the density of states for the pure system N(E) by the averaged 
one of the doped material N(E). Evidently, assuming small fluctuations in the the pairing 
potential (Eq. |2S[) one gets a critical temperature T c weakly depending on disorder. This 
is the content of the Anderson theorenJH which rests on the assumption that pairing 
potential does not fluctuate in space A» ~ A for all i. However, it should be noted that 
in short coherence length limit the situation can be opposite. In that case allowing the 
spatial fluctuations of the pairing amplitude Aj ^ Aj induced by the site energy disorder 
Ei 7^ Ej or introduced by the randomly distributed attractive centers (Uu = for some 
lattice sites i) the Anderson theorem breaks down even for conventional, on site, s-wave 
superconductorsB0. 

For a constant pairing parameter as in Eq. ( p8|) the generic consequence of disorder in 
the system with on site attraction is the smearing of the structure in the averaged density of 
states N(E) (Eq. [£|). To illustrate the consequences of it in our model we have calculated 
N(E) using the standard CPA procedure@S0 for c = 0.5, Ea = — f , £b 



_ 8 

2 , 2 - 
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The results for various values of the scattering strength 5 leading to the smearing of the 
Van Hove singularities in the averaged density of states N(E), are shown in Fig. 5a. Van 
Hove singularities are still present here for a relatively weak disorder strength 8 < It while 
for stronger one 5 = 2t one can notice additional splitting of singularities caused by the 
model of disorder. This is the so called the split band regimeBEl. Namely the two peaks 
are the remnant of the Van Hove singularities of the two, A and B, pure metals. 

Let us now examine a disordered system with the inter-site attraction C/y. Here we 
assume that inter-site pairing parameter Ay can be substituted by its average 
Thus, in the case of a singlet paring (extended s- or d- wave), the on site impurity potential 
V a can be expressed as 



V ( 



e a 
-e, 



a = A or B 



(35) 



while the coherent potential matrix has the form: 



£(£) 



E u (£) 
£ 22 (£) 



(36) 



Naturally, the averaged Green function, G(i,i; E) can be expressed as follows: 



G(y;i?) = lEG(M)4E 



E-e k -E n (E) A fc 

A; E + t k - Y, 22 {E) 



(37) 



In case of triplet pairing (p-wave) instead Eqs. (|35"l - j37j ) the following notation should be 
introducedii'@ 



e a 

e a 

-e a 

-e a 



a = A or B 



(38) 



for the impurity potential and the coherent potential 
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£(£) 



En (£7) 

En(-B) 

E 22 (£) 

E 22 (£) 

respectively. Again the averaged Green function is given by 



(39) 



N 



-i 



(40) 



(E-e k -E n (E))l A fc 
k ' k A; (£ + e fc -E 22 (£))l 

while the conditionally averaged Green function at the impurity site (a = A or 5) has the 
following form: 



G a {i, i; E) = G(i, z; £)(1 - [V a - X a (E)]G(i, z; E))~ x . 



(41) 



The averaged paring parameters can be written as in (Eqs. 113): 



K = Ao7 fc , 
Afc = A r/ fc , 

A^ = AQsin(fc ;t a) + A^ sin(k y a) 



(42) 



For anisotropic s— , d- and p-wave pairing symmetries the gap equations (Eqs. |9| and [10D 
take the form: 



Afc = ^ E r &E lmG 12 (k; E) ' 



iV g 7T J-oo 

for singlets (extended s- and d- wave) and 



+ 1 



(43) 



(44) 



for triplets (p-wave) cases. 

From Eqs. ( |3~TD and ([!(]) it follows that off diagonal elements of the Green function are 



G u (k;E) + G 22 (k;E) ^ 
G ^ E > ~ 2£-S 11 (£)-E 22 (£) A ° Cfc ' 



(45) 
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for singlets and 

G 12 {k;E) = 2E - Z n (E) - Z 2 2{E) ^ Sm Sm ^ ' ^ 

for triplets, where ( k = 7 fc or r] k . 

Interestingly, the linearized gap equation can be written as in a way similar to that for 
the clean system (Eq. 

1 = Hr««| ta ^-. (47) 
7T y_oa 2 2£-Tr£(£) v ; 

where the imaginary parts of G s,d,P (E) define the projected densities of states for a disordered 
system N S (E), Nd(E) and N P (E) discussed in Appendix A. Moreover, as in case of clean 
system (Eqs. p0|) : 

N s , d , P (E) = ~lmG sAs (E) = --^ECym 1 (48) 

k ~ ^n\E) — e k + fi 

where (' k = (7fc) 2 /4, (r/ fc ) 2 /2 or 2 (sin (A^ a)) 2 depending on the symmetry of solution: ex- 
tended s-, d- and p-wave. Figures 5b-d show the corresponding projected densities of 
states of a disordered system (Eq. g|): N S (E), N d (E) and N P (E). Like for N(E) (Fig. 5a) 
for weak disorder, Van Hove singularities survive. Thus, it is clear that the smearing of the 
Van Hove singularity in Figs. 6a-d implies a weakening of the Van Hove enhancement of the 
critical temperature T c . 

For our simple model of disorder of binary alloy A c Bi_ c with c=0.5 the coherent potential 
Yiii(E), satisfies the following CPA equation^!: 

E U (E) = (±8 - Z u (E))G u (z, i; E)(^S + E U (E)). (49) 

In Fig. 6a,b we also show the corresponding self energy £(£") which in one site CPA 
depends only on energy E but not on the wave vector k. One can see that for a relatively 
weak disorder strength S the maximum of |ImSu(£J)| is exactly at Van Hove singularity as 



in the Born approximationlllillS 



r 2 

ImS(E) « -^N(E). (50) 
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The above result is valid also below the critical temperature T c because for off diagonal 
pairing A u = oEHH and (Eq. |9|) is still valid. 

In case with a larger disorder strength 5 the maximum of |ImSn(i?)| is located in some 
other energy according to the model of disorder we use (Eqs. [I], |35| - f4TD . Simultaneously 
the real part of self energy ReSn(i?) is changing with energy E renormalizing the chemical 
potential /i (Eq. |48|) . 



Let us investigate the additional effect of disorder visible in Eq. (f47|) . Comparing this 
equation with the clean system one (Eq. [19|) one can notice the difference in the denominator, 
where in case of disordered system there is an additional strong scattering term Tr S(i?). 

To examine it further let us rewrite the gap equation in terms of Matsubara frequencies 



(3 C V tuJ n - En(zo;„)' 

for s- d- and p-wave symmetry. 

Now, let us approximate the normal state self energy Sn(zc(j„) and projected density of 
states N a (iu> n ) by: 

£(m; n ) « -«|E |sgn(a; n ) (52) 
N a (uo n ) = N a {wj n - S(zo; n )) w N(u n + 2|S |sgn(^„)) (53) 
= —ImGx^i, i; iu n + «|S |sgn(u; n )) 

7T 

In Fig. 7a we have plotted the densities N d {iui) N p {iuj) N s {iui) versus imaginary part of 
energy iuj for Fermi energy Ef chosen at Van Hove singularity (Ef = E v = 0). 

Note that in that region we can roughly approximate the corresponding projected den- 
sities by a simple formula 

N a (tu n ) = — ^V, (54) 

where a a and b a are constants depending on the pairing symmetry a = s,d,p. Introducing 
Eq. O and Eqs. (HH) into Eq. (0) we get: 
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1 = \U\a*T e £ x — (55) 

Now we have to perform the summation over u n for clean Eo = and disordered cases 
Sq 7^ 0. After some algebra (Appendix C) we get the approximate pair-breaking formula: 



*G) 



I b n \ , ( 1 S \ / 1 S 6 n 



n^^t ^u+^-n^^ ■ (56) 



Note that in the limit the constant density of states N a {iuj) = const. (Eq. [54]), b a — > oo, 
we get automatically the the standard Abrikosov-Gorkov formula! with a characteristic 
pair-breaking parameter p c = |Im£o|/(27rT c ): 



1» ( ^ 

T, 



cO 



*G)-*G + *) (57) 



In Eqs. and |57| ) T c0 denotes the critical temperature in a clean system, while T c is the 
critical temperature in dirty one. 

In Fig. 7b we plot T c /T c0 versus S /(27rT c0 ) for few values of b a . Here the chemical 
potential /i was fixed at the saddle point Van Hove singularity E v \. Interestingly, the slop of 
the curve (Fig. 7b) is increasing with decreasing b a indicating that in presence of Van Hove 
singularities for a disordered system we get weaker decreasing of T c than for a flat density of 
states. Thus, for d-wave pairing b a = 0.12 and for p-wave b a = 0.45 superconducting phase 
is more stable than in case of extended s-wave pairing, where b a — > oo (Fig. 7a). This is 
the main result obtained in this section. In spite of very rough approximation used here the 
results show that the Van Hove singularity influences the pair-breaking effect. Namely it 
makes superconductivity more robust increasing the critical strength of disorder E needed 
to break the Cooper pairs. 



IV. CRITICAL TEMPERATURE FOR DISORDERED SUPERCONDUCTORS 



Let us now turn to the case where both superconductivity and disorder are presentt^il^ 
and calculate the critical temperature T c self-consistently. Although the full CPA program 
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can be implemented for the problem defined by Eqs. (|9|-|r0| and [B.l|4B.4|) Ellzl and the speci 



fication of the site energy ensemble, it is convenient to make the approximation, valid when 
the coherence length £o is much larger then the lattice spacing, that the pairing potential 
Ay does not fluctuate very much and replace it in equations (| and ||) by its average value 
A„§. For conventional isotopic .-wave pairing the gap equation at T c takes the simple form 
(Eq. [34]). Thus, the critical temperature T c at the optimal doping should be only slightly 
reduced by disorder due to smearing of the density of states N(E). The results of numerical 
calculations for four different values of a disorder strength 5 in case of the 'on site' attraction 
is presented in Fig. 4c. Clearly, in this case the critical temperature T c is sightly reduced 
by the effect of the density of states. 

On the other hand for the 'off diagonal' attraction case the linearized gap equation, in 
presence of disorder, is given by the formula (Eq. ^7|). Solving it for extended s- d- and p- 
wave pairing symmetries we get the critical temperature Tc versus band filling. The results 
for various 5 are presented in Fig. 4a (extended s- and d-wave) and 4b (p-wave). Here, 
we observe significant degradation of T c in all three anisotropic pairing cases. Moreover, 
in some regions of electron concentration, disappearing of a supercondacting phase can be 
noticed for relatively weak disorder (5 < 0.6t). Especially this happens to the s-wave case 
with high electron concentration n — > 2 and n m 0.8 as well as the <i-wave case for n ~ 1.2 
(Fig. 4a). The dramatic degradation of extended s-wave superconductivity for n > 0.3 can 
be also seen for stronger disorder (5 = It). In that region of electron concentration n, due 
to the binary alloy model of disorder A C B!_ C and c = 0.5, the position of energy E for the 
maximum of a pair-breaking term |ImEn(.E)| « So (Fig. 6) coincides with the chemical 
potential /i making the pair-breaking mechanism very efficient. Similar behaviour can be 
seen for p-wave superconductor (Fig. 4b). Clearly, for large enough |ImSn(i?)| (Fig. 4a-b) 
the superconductivity is destroyed by the pair-breaking effect shrinking the region of n for 
T c > 0. In the same time the corresponding projected densities are not strongly effected 
(Fig. 5). The results for anisotropic paring in Figs. 4a-b are contrasting with Fig. 4c where 
we plotted the results for the on site solution with U = Ua. Here, (Fig. 4c) the region of 
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band filling n, where a superconducting solution exists, does not change with disorder at all. 

V. CONCLUSIONS AND REMARKS 

We have analyzed the effect of disorder on disordered Hubbard model with local and 
non-local, nearest neighbour, interactions as well as nearest and next neighbour electron 
hopping terms. We have got numerical and analytical results confirming previous works on 
the similar model with a simple band energyElS e& = — 2t(cos k x a+ cos k y a). Such dispersion 
relation, introduces the electron-hole symmetry for half filled band n — 1 and locates the 
saddle point Van Hove singularity exactly in the center of the band (Figs. la,b). Including 
the additional hopping t', we break the electron hole symmetry in the densities of states, 
resulting in shifting the central Van Hove saddle point singularity to the bottom of the band. 
In the same time another Van Hove singularity, located at the bottom band edge, appears 
to be important. Thus the effect of Fermi surface distortion make place both of singularities 
very close to each other For some of band filling values n « 0.4 both of singularities play 
important roles. It is clear after analyzing the appropriate projected densities of states 
N a . Moreover various pairing symmetries choose different Van Hove singularities. D-w&ve 
symmetry is favored in case of the regions of band fillings n with chemical potential \x near the 
saddle point singularity E v while extended s- wave pairing symmetry is more likely as far as 
the bottom edge singularity is concerned. Interestingly, on account of the non-symmetricity 
of electron density of states the Van Hove scenario is fulfilled only approximately. Here we 
observe a small interaction dependent, shift of the optimal doping electron concentration 
n m n op towards the center of the band. As a result of the above we conclude that Van Hove 
singularity is important for all discussed symmetries of order parameter. Interestingly, in 
presence of an additional electron hopping to the next nieghbour lattice site t', the shift of the 
saddle point Van Hove singularity in the density of states explains that the superconductor 
with an extended s-wave symmetry is rather of electron type while with ci-wave symmetry 
case is of hole nature. Note also that, the positions of the maximum value of N S (E), Nd(E), 
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N p (E) (Fig. 5) are not affected by small disorder. In spite of small density of states effects 
leading to smearing peaks in corresponding densities: N S (E), Nd(E) and N P (E) (Eq. |48|, 
Fig. 5), the critical temperature T c , plotted in Figs. 4a-b, is degradated strongly with 
disorder. This is due to the pair-breaking term Hu(iu) n ) (Eq. |51]). In the negative U, 'on- 
site', interaction there is a quite different situation. Here disorder causes only small decrease 
of T c (Fig. 4c) via a density of states effect (Fig. 5a). So, the most interesting effect arises 
from Eqs. (|51|j55| ), where En (£7) acts as a pair breaker. 

Concluding our results we would like to stress that the Van Hove singularity does not 
make the decrease in T c more pronounced as one can expect (because of singularities are 
present in the self energy ImSn(E) ~ N(E), Fig. 6). In fact singularities make it weaker. 
We have analyzed this effect very carefully in Sec. Ill finding an approximate pair-breaking 
formula (Eq. [56]). Here, the Van Hove singularity influenced the pair-breaking Abrikosov- 
Gorkov curve, changing its slope of the T c versus E . Similar effect has been observed in the 
experimental results for Zn-doped LSCOll. Alternatively, this effect can be also explained 



assuming anisotropic scattering potentialsBQ. 

Finally, we observed the dependence of T c on band filling n. Our results for d-wave 
superconductor show that the Van Hove scenario is valid even in presence of weak disorder 
(Figs. 4). Similar experimental results were obtained by measuring T c , in various cuprate 
compounds, as a function of hole concentration where Cu were substituted by Zn 0H. The 
concentration of Zn was there the measure of disorder. 

It should be however noted that high T c cuprates are strongly correlated electron systems 
and the mean field approach basing on the effective intersite attraction Uy, presented here 
has a limited applicability!!. Strictly speaking a more realistic model should posses, beside 
an intersite attraction, a strong Coulomb repulsion term. Although the approximations we 
used in the present paper aimed to explain the effect of the Van Hove singularity in a presence 
of disorder in a weak coupling regime, the preliminary calculations using simultaneously 
slave boson technique and CPA0 seams to support the general arguments of the Van Hove 
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singularities significance for superconducting cuprates conjectured here. 



A. APPENDIX 

In this Appendix we apply the recursion method to calculate appropriate densities of 
states. Let us investigate the projected densities if states N S (E), Nd(E), N P (E) and the 
corresponding Green functions: G S (E), G d (E) and 

N a (E) = --lmG s (E) = -i-Ey -ImGni^E), 

7T 1\ T 4 7T 

N d (E) = --lmG d (E) = — E T ~ lmG ^ E )> (A-l) 

7T iV ^ 4 7T 

iV p (£) = --ImG p (£) = -^2(siny 2 -lmG n (k,E), 

where 7 fc and rj k were defined in Eq. 11. 
Noting trigonometric identities: 



t| i 

4 2 



— = -(cos 2k x + cos 2k y ) + 1 + 2 cos k x cos k yi 



rf 1 

— = - (cos 2k x + cos 2k y ) + 1 — 2 cos k x cos k yi (A. 2) 

2(sinA; i .) 2 = 1 — cos2/c x , 

the Green functions G S (E), G d (E) and G P (E) can be easily found as a combination of 
diagonal and off diagonal Green functions Gn(i + 5i, i + 5j, i?) = Gsi,sj{E), see the notation 
in the schematic picture (Fig. 8): 

G S (E) = G 00 (E) + G 20 (E) + 2G U (E), 

G d (E) = G 00 (E) + G 20 (E) - 2G U (E), (A.3) 
G P (E) = G 00 (E)-G 20 (E). 

Function Gqo(E), G 2 q(E) and Gn(E) have been calculated using the recursion methodic. 

The above procedure can be also used to calculate the average of projected densities 
N S (E), Nd(E) and iV p (.E)as well as averaged Green functions G S (E), G d (E) and G P (E) can 
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be calculated via substitution E by E — T,u(E), where local £n(J5) = Sn(i, i, E) should be 
found self-consistently according to CPA conditions (Eqs. |37| - |43j ) : 

G sAp {E) = G sAp (E - E n (E)) 

N sAp (E) = N sAp (E - £ u (£))- (A.4) 



B. APPENDIX 



In this Appendix we apply CPA equations (Eqs. |18|-^1|) to the disordered Hubbard model 
with the 'on site attraction' Eq. [I] (U^ = Uu5ij) and discuss the Anderson theoremS'El. The 
averaged Green function, G(i, i; E) can be expressed as follows: 



£-e fe -£ u (£) -S 12 (E) 
-E 21 (E) E + e k -X 22 {E) 

while conditionally averaged Green function has the form 



G^;E) = ±Y.G(k;E) = ±Y: 

k k 



(B.l) 



G a (i, z; E) = G(i, z; E)(l - [V a - E a (£7)]G(i, z; S))" 1 . 

The disordered potential in Eq. 18 has the following form&fl: 

e a -A c 
-A* — F 



(B.2) 



(B.3) 



where e a = Ea or e b corresponds to different site energies of the lattice site while A a = A a 
or As relates to the different paring potential in an alloy A C B\_ C . 
Clearly, the coherent potential can be written as 



£ n (£) E 12 {E) 
Z 21 (E) £ 22 (£) 



(B.4) 



Function G\ 2 (i,i; E) can be obtained from Eq. [B.l 



3 = G^EHG^E)^ = /^(m;^ + ^,(m;JB) \ 



2£-E 11 (£)-E 22 (£) 



2£-E u (£)-E 22 (£) 
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While from Eq. p.4] we have the following relations: 



TrG«(*, i; E) = DetG a (i, i; E) { ^^^fl) ~ ^ ll{E) ~ ^ 22{E) ) (B ' 6) 
GUh »; E) = DetG<*(i, i- E) (^^^y - A« - X 12 (E)\ (B.7) 



substitution of Eq. |B.6| and |B.7| into the right and left hand sides of Eq. B.5| respectively 



an equation on E^-E 1 ) and A a : 

2E S 12 (E) (A a DetG a (i } i; E)) 

2E-Z n (E)-Z 22 (E) ~ (T)etG(i,i;E)) 

Using |B.6| and factorizing of above we get: 

2EZ 12 {E) cTiG A {i,i;E)A A + (1 - c)TrG B E)A B 



(B.8) 



2E - S n (E) - S 22 (E) TrG(i, i; E) 

In the limit of small paring potential fluctuations it can be used as a constant A, — > A 

const, to satisfy the Anderson theorem. Then 

2E E 12 (E) 



B.9) 



2E-E U (E)-Z 22 (E) 
and renormalized quantities of ui and A reads as: 

2E = 2E~ TrS(E) 



A (B.10) 



^= W = ™>S (B.ll) 



This lead to the same renormalization in the pair potential A and energy E: 

A(E) _ E 
A ~ E 



(B.12) 



A = — — f°° d E MG^ + G^A 1 

7T J-oo 2E e? E + 1 v ; 

U r iE Im5 11 (^)A tajih 



7T J-oo 2£ \ 2 / 

and the linearized gap equation 



1 = 17 / F(E)d£ tanh(/ ^ (B.14) 



where N(E) denotes normal state DOS: 



JV(£) = — ImGii(^). (B.15) 

7T 
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C. APPENDIX 



In this Appendix we apply CPA to anisotropic superconductor and evaluate the approx- 
imate formula of the pair-breaking effect. Starting from the clean system we assume that 
the linearized gap equation can be written (Eqs. [53|-^7]) by: 



1 = \U\a^T c0 £ ^- (C.l) 

\U\a a 2irT c0 ^ / 1 

~ — i — E — 



Ja w„>0 



where uj c n is a cut-off in the summation ( |C.1| ) As io c n is very large, it yields 



On the other hand for disordered system we have: 



1 = \U\a^T c J2 * v | x — i (C.3) 



Similarly (to Eqs. |C.1| , |C.2| ) it leads to 



" U + Pc ; + Pe + 7rlr), (C4) 



|E/|a« V2 r V V 2 27rT - 
and finally comparing Eqs. (|C.2| and |C.4|) we get 



where p c is a pair-breaking parameter 



27rT t 

and T c0 is the critical temperature for a clean superconductor 



Pc = S (C.6) 
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FIGURES 



LU 




FIG. 1. (a) The distortion of the Fermi surface by the next neighbour electron hopping t' . The 
full line denotes Fermi surface for if = while the dashed line if = 0.451 Fermi energy is at Van 
Hove singularity Ep = E v . (b) The electron density of states N{E) for the 2D lattice, where 
the full line corresponds to t' = and the dashed line to t' = 0.45i. Arrows denote Van Hove 
singularities E v , E v \ and E v i, respectively (Ep = 0). 
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FIG. 2. The band structure (a) and Fermi surfaces (a,b) for the one band electron structure 
with next nearest neighbour hopping: e& = — 2i(cos k x + cos k y ) + At' cos k x cos k y with t' = 0A5t, 
and three different band fillings: n = 0.584 (1), n = 0.467 (2), n = 0.388 (3). 
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FIG. 3. (a) The electron density of states N(E) (full line) and projected densities (dotted 
lines) for extended s-wave N S (E) and ci-wave type Nd{E) (a). N(E) (full line) and p-wave type 
N p (E) (dotted line) for a normal pure system, the chemical potential \i = 0. 
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FIG. 4. (a) The critical temperature T c versus band filling n for extended s-wave and d-wave 
pairing for a pure '1' and alloyed disordered system '2'-'4' (depending on 5); and U = — 1M. Full 
lines corresponds to d-wave solution while dashed extended s-wave ones, (b) T c versus band filling 
n for p-wave solution, (c) For comparison T c (n) for on site s-wave solution. 
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FIG. 5. Normal state density of states N(E) (a) and projected densities: extended s-wave 
N S (E) (b), d-wave type Nd(E) (c) and p-wave N p (E) (d), respectively, for different values of 
disorder strength S/t = 0.0, 0.1, 0.2, 0.3. Arrows show the directions of S change. Here, the 
chemical potential fi = 0. 
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FIG. 6. Normal state parts of self energy S(-E): imaginary (a) and real (b) for different values 



of 5. 
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FIG. 7. (a) The imaginary part of Green functions: — ImG^, — ImG p and —lmG s versus 
imaginary energy iu. Fitted lines a a /(b a + uj) are plotted for d-wave pairing where b a = 0.12 
and p-wave pairing b a = 0.12 (a a = 2/3 for both curves), (b) The critical temperature T c 
versus pair breaking parameter |Eo| in the limit of weak disorder (both T c and |So| are normalized 
to T c q of clean system) for few values of b a . The limit b a — > oo corresponds to the standard 
Abrikosov-Gorkov formula (Eqs. 58-59). 



39 



o o o o o 



o 




o 



o 



X 



FIG. 8. Schematic picture of diagonal and off-diagonal Green functions G Q / g = G(ax + fly). 
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